Neutralino Relic Density in the CPVMSSM at the ILC 
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We discuss ILC measurements for a specific MSSM scenario with CP phases, where the 
lightest neutralino, a candidate for dark matter, annihilates through t-channel exchange 
of light staus. These prospective ILC measurements are used to fit the underlying model 
parameters. A collider prediction of the relic density of the neutralino from this fit gives 
0.116 < flh 2 < 0.19 at 95% CL. 



1 Introduction 

One of the prime motivations for a high-luminosity e + e~ linear collider is the possibility to 
do precision measurements of new particles beyond the Standard Model and in particular 
of the dark matter (DM) candidate. A precise determination of the properties of the DM in 
the laboratory could be used to make a "collider prediction" of its relic abundance, which 
can be tested against cosmological models. For such a collider prediction to be of interest, 
it must be at least as precise as the value obtained from cosmological observations, about 
10% for WMAP+SDSS [2] and a few percent at the PLANCK satellite [3]. 

The possibility to make collider predictions of the cross sections for annihilation of dark 
matter candidates has been examined within specific supersymmetric scenarios. For ex- 
ample, in the general MSSM, it was shown that in a favourable scenario the LHC could 
match roughly the WMAP+SDSS precision [4j while the ILC could achieve much better 
precision [5]. These conclusions, however, depend very strongly on the scenario considered; 
many remain challenging even at the ILC [5]. 

These studies assumed that CP is conserved, although CP-violating (CPV) phases are 
generic in the MSSM. CPV MSSM phases can have an important effect on the neutralino 
annihilation cross sections [5] and could lead to variations in flh 2 solely from modifications in 
the couplings of up to an order of magnitude [7] . Therefore the determination of the relevant 
couplings (including phases) can be as important for the prediction of flh 2 as measurements 
of masses. In a particular scenario of the CPV MSSM, the precision to which the underlying 
MSSM parameters could be determined from measurements at the ILC was investigated 
in [8]. The resulting presicion on the neutralino relic density that could be inferred from 
this was examined. The main reults are summarised below. 
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2 The stau-bulk scenario of the CPVMSSM 



The scenario we investigate belongs to the "stau bulk" region of [7] , which appears for light 
staus and large phase of Mi, the gaugino mass. The neutralino LSP is dominantly bino and 
annihilates predominantly into tau pairs, via t-channel exchange of both t\ and f 2 . The 
annihilation cross section is sensitive to the stau mixing . Although the f\ is very light, the 
scenario does not rely on coannihilation (the mass difference with the LSP is too large). 
As a benchmark point we choose the following input parameters at the electroweak scale: 

Mi = 80.47 GeV, M 2 = 170.35 GeV, M 3 = 700 GeV, fa = n, fx = 600 GcV , . 
M La = 138.7 GcV, M &3 = 135.2 GcV, A T = 60 GcV, fa = 0, tan/3 = 10. 

Other sfcrmion masses, A t .b, and m H + are set to 1 TeV and </> M = 0. This way EDM 
constraints are avoided when varying fa and fa. The particles accessible at ILC500 are x?, 
xi,Xi, n, f 2 and v r with masses respectively of (80.7,164.9,164.9,100.9,177.2,123.1) GcV 
as computed with CPsuperH [5]. Despite having a light spectrum this scenario is quite 
challenging for colliders. Production of sparticles at ILC500 only lead to r's plus E^ lss and 
similarly for the cascade decays of squarks at LHC. 

The relic density of the Xi computed with micrOMEGAs2 . 2 10J is flh 2 = 0.130. The 
precision with which fill 2 can be inferred from ILC measurements depends not only on the 
accuracy of the sparticle spectroscopy but also on the determination of all parameters of the 
neutralino sector (Mi, M 2 , /J,, tan/3, fa) and of the stau sector {Mf L , Mf R , A T , fa). The 
dependence on fa, and to a much lesser extent fa, originates from the Xi^i i T couplings. In 
addition, particles which are too heavy to be produced at ILC could have some influence. 
We assume here that the mass scale of the squarks and gluino is known from LHC and that 
the ILC excludes selectrons/smuons and heavy Higgs states up the the kinematic limit. This 
leaves an overall uncertainty from the unknown part of the spectrum of S^l/Sl ~ 7%. 

3 ILC measurements 

At the ILC, sparticle masses can be mea- 
sured either with threshold scans or with 
the endpoint method. The challenge in 
this scenario is to disentangle the various 
sources that lead to the the tt£"" ss topol- 
ogy. Scanning downwards in energy in steps 
of 10 GeV while using different beam polari- 
sations allows to detect thresholds: the <jrl 
polarisation mode for f$ ~fj pair production 
and the <7lr mode for x + X~ production. In 
addition the mixing angle in the stau sector, 
cos2#f , can be determined from measurements of the polarized cross section while the mea- 
surement of the r polarisation V T in the decay f\ — ► tx± gives additional information on the 
stau and neutralino mixings |11| . It is in particular useful to constrain the gaugino-higgsino 
composition of the LSP. It was shown in [8] that masses, polarisation V r and mixing 9f can 
be accurately determined even with moderate integrated luminosity, see Table [TJ This anal- 
ysis used simulations based on Simdet 4.02 [12], beam polarisations ('P e -,7 : ' e +) = (0.8,0.6) 
and includes QED radiation, beamstrahlung and SM background to W pairs. 



channel 


observables 


ri + ff 


ijifi = 100.92 ± 0.40 GeV 




m^o = 80.67 ±0.35 GeV 




cos2 f = -0.065 ±0.028 




V T = 0.64 ± 0.035 




m f2 = 176.9 ±9.1 GeV 


XiXi 


m ± = 164.88 ±0.015 GeV 



Table 1: Achievable precisions at ILC. 
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4 DM properties: fit to ILC observables 



We performed a fit to the ILC observables listed in Table Q] and used a Markov Chain Monte 
Carlo (MCMC) method to probe the 8-dimensional parameter space. The free parameters 
are Mi, \x, tan/3, M^, M^, A T , (f>\, (f> T , while Mi is computed from fj. and tan/3 to match 
the extremely well measured value of m^± . The total tt SUSY cross section with polarized 

beams was also included in the fit. Moreover, we have a TT = 3.220 ± 0.046 pb at y/s = 
400 GcV including the systematic and statistical uncertainties for 10 fb _1 . 

The very precise determination of m^o and 
m-+ constrain Mi to SMi = -0.9, -^.l'GeV 
and indicates a correlation between §i and Mi 
although the phase 4>i can take any value. The 
measurement of the stau masses and their mix- 
ing angle constrain the stau soft masses to 
about 10 GeV. The trilinear coupling A T and 
its phase 4> T are on the other hand basically 
undetermined. This is because the stau mixing 
which is proportional to ( A T — /i tan /3) is domi- 
nated by the term fj, tan (3. This also means that 
fi and tan/3 are individually poorly determined, 
see Fig. [TJ 

Within the allowed parameter space, the 
relic density of dark matter is predicted to be 
0.116 < nh 2 < 0.19 at 2a. The largest values 
of Vlh? are found for <fii « and small \i, where 

the LSP has an increased higgsino component, see Fig. [5] An improvement on the limit on 
the neutralino-proton cross section, which is largest when /i is small, would cut on these 
scenarios thus somewhat reducing the range for ilh 2 . 

A better determination of tan/3 and \i would be needed to significantly improve the 
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Figure 1: 2a (green) and la (black) al- 
lowed region in the tan/3 vs /i plane. 
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Figure 2: Predictions for flh 2 vs a) fj, and b) 4>i and the WMAP upper limit (line). 
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uncertainty in Vlti 2 . For both, a multi-TeV e + e~ linear collider offers the best prospects: 
The parameter /i could be determined from a measurement of the heavy higgsino states. 
A measurement of the stop masses and mixing angle at LHC and a multi-TeV LC would 
completely determine the stop sector thus reducing the large parametric uncertainty in rrih- 
The precise determination of mt at the ILC could then be used to constrain tan/3. 

CP-odd observables such as electric dipole moments or T-odd asymmetries could give 
a very clear signal of CP violation. However, it was shown in [5] that a T-odd asymmetry 
based on triple-products in the production and decay of neutralinos [13] while sensitive to 
tfii , would not help constrain flh 2 . The reason is a two-fold ambiguity in <f>i . 

5 Conclusions 

A precise determination of SUSY particle properties in the "stau-bulk" scenario of the 
CPVMSSM is quite challenging at the ILC. Although some of the underlying Lagrangian 
parameters can be extracted with very good precision, the large uncertainties in tan/3 
and the phase <pi induce a large uncertainty in the prediction of the neutralino relic density. 
To reduce this uncertainty to the percent-level would require precision measurements of the 
heavy higgsino-like neutralinos and charginos at TcV energies. 
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